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Abstract. Gerbillus campestris is a widely distributed small rodent that lives in various habitats in North Africa and can be 
a potential agricultural pest. We conducted a biogeographic study of this species with an integrative approach using mor- 
phometric data from body and craniomandibular distances and molecular genetic data from the cytochrome b gene of the 
mitochondrial DNA. We collected 96 individuals in six localities from central, northern and eastern Morocco. Data from 
18 morphological characters were used in multivariate statistical analyses and molecular data were analysed using maxi- 
mum likelihood and median-joining networks. Our analyses confirmed the high morphological variability in this species 
and allowed to discriminate four groups containing the studied populations. We found that a few craniomandibular measu- 
rements had the highest contribution in the differences between populations, and that this variability reflected a spatial and 
environmental differentiation. In the genetic analyses, we placed our six populations and six individuals from four other 
countries in nine previously identified phylogenetic lineages in this species, and we also added a tenth lineage. Limited 
gene flow, isolation by distance and biogeographic barriers were further explored to explain this genetic structuration. We 
also jointly examined morphometric and genetic variability and found that the morphological groups were congruent with 
the genetic lineages and the geographic distribution. A better knowledge of the phenotypic plasticity and genetic diversity 


of this gerbil can be used to comprehend the micro-evolutionary processes in other small mammals in North Africa. 
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INTRODUCTION 


The Gerbillinae subfamily and the genus Gerbillus are 
the most diverse group of rodents in North Africa. Like 
other small rodents, the species of this genus represent 
suitable models for the construction of biogeographic his- 
tory thanks to their brief generation time and prolific re- 
production rates; their rather limited dispersion capacities 
and strong affiliation to arid and/or open habitats (Avise 
2000; Nicolas et al. 2009). A few recent molecular studies 
have attempted to review the systematics of the various 
species of Gerbillus inhabiting North Africa and Moroc- 
co, or to better understand their biogeography (Abiadh 
et al. 2010; Ndiaye et al. 2012, 2013a, 2016; Nicolas 
et al. 2014; Bouarakia et al. 2018). Among this genus, the 
most widespread species in Morocco is the North African 
Gerbil, Gerbillus campestris Loche, 1867, which is also 
considered as an agricultural pest (Giban & Haltebourg 
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1965; Ouzaouit 1980; Zyadi & Benazzou 1992). It is a 
long tailed, middle sized gerbil, with naked hind feet, that 
lives throughout the Mediterranean part of North Africa 
from Egypt, west of the Nile River and Delta, to Morocco 
in the West. It is also present but less widely distributed in 
the Sahara Desert and Sahel region, in Sudan, Niger and 
Mali (Happold 2013; Granjon et al. 2016; Denys et al. 
2017). The species dwells in different habitats, from sub- 
humid to arid and desertic regions, except high moun- 
tains, forests and sandy desert areas. It occupies steppes, 
arable lands, rocky habitats and oases; and digs its bur- 
row in sandy or clay soils (Petter & Saint Girons 1965, 
Happold 1967; Aulagnier & Thevenot 1986; Kowalski & 
Rzebik-Kowalska 1991; Dobigny et al. 2002; Musser & 
Carleton 2005; Aulagnier et al. 2017). 

Previous studies have investigated the intraspecific 
variability of Gerbillus campestris. Petter & Saint-Gi- 
rons (1965) suggested the subdivision of the species in 
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Morocco into to two subspecies based on coloration cri- 
teria. In Tunisia, due to geographic variation in size, the 
existence of two subspecies was suggested (Ranck 1968; 
Jordan et al. 1974). By combining the variation of skull 
size and characters, external size and pelage colour, two 
to five subspecies were identified in Egypt, Libya and Su- 
dan (Setzer 1958; Ranck 1968; Osborn & Helmy 1980). 
Benazzou & Zyadi (1990) and Baala (1995) displayed 
the presence of a geographic biometric variation between 
five populations in Morocco, while Baala (1995) showed 
the absence of biochemical polymorphism between these 
five populations. On a chromosomal level, the most com- 
mon diploid number of chromosomes for G. campestris 
is 2n = 56. This karyotype is present in Algeria, Tunisia, 
Egypt, Mali and Niger (Matthey 1952; Wassif et al. 1969; 
Jordan et al. 1974; Dobigny et al. 2001, 2002). However, 
chromosomal polymorphism was documented in Moroc- 
co with 2n = 56 in Taroudant locality and 2n = 57, 58 in 
Goulimine locality (Lay et al. 1975). More recently, Nico- 
las et al. (2014) presented a complex phylogeographic 
pattern in Morocco and North Africa using cytochrome 
b sequences. The wide distribution of G. campestris, the 
multiplicity of habitats it occupies, the rich geographic 
diversity of the region and the extreme environmental 
changes that occurred in North Africa during the Pleisto- 
cene, parallel this important intraspecific diversity. 

To complete the phylogeographic structuration uncov- 
ered by Nicolas et al. (2014), we increased sample size 
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in three localities and added three new ones in Morocco, 
to perform a study of genetic variability based on the cy- 
tochrome b gene of the mitochondrial DNA. The addi- 
tional sampling from north-western Morocco helped us 
to explore the existence of barriers to gene flow in this 
part of Morocco, namely the Sebou River. In addition, 
the new sampling from north-eastern Morocco on both 
the eastern and western sides of the Moulouya River was 
useful to investigate the role of the Moulouya River and/ 
or its arid valley as a major geographic barrier and to ver- 
ify the suggested east-west genetic differentiation around 
this river (Nicolas et al. 2014; Lalis et al. 2016; Beddek 
et al. 2018). We also included samples from GenBank 
from four new localities (Mali, Tunisia, Libya and Egypt) 
to better identify past gene flow in other parts of the dis- 
tribution range of the species. 

Furthermore, we used body and craniomandibular dis- 
tance-based measurements of genetically typed individu- 
als from six localities from central, northern and eastern 
Morocco in a comparative study of morphometric vari- 
ability, to explore the influence of spatial, geographic and 
environmental factors on phenotypic diversity. Finally, 
we jointly examined the morphological and genetic inter- 
population variability to verify their congruence in defin- 
ing the biogeographic history of the species in Morocco. 


Fig. 1. Map displaying the geographic range (grey area, as mentioned by Granjon et al. 2016) of Gerbillus campestris and the dis- 
tribution of the ten mitochondrial lineages. Dotted lines surround the geographic groupings of five lineages (1, 2, 3, 4 and 6). Based 
on Nicolas et al. (2014), only the localities of Mali and Niger were included in the geographic grouping of lineage 2. Localities from 
the present work (CHR, ANT, KTR, MZ, GCF, FRT) and from bibliography (numbered from 1 to 28) are indicated (see Appendix 
I). Localities 14, 16, 20, 23 and 24 shelter haplotypes from both lineages 1 and 3. Moulouya, Sebou, Niger and Nile Rivers are also 


shown in the map. 
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MATERIALAND METHODS 


Sampling. We collected 96 specimens of Gerbillus 
campestris between 2011 and 2015 in six localities in 
Morocco from various regions: 10 specimens in the 
Tingitane peninsula in Chrouda locality (named hereafter 
CHR); 13 in central Morocco in Aounate locality (ANT); 
24 in Kenitra (KTR) and 25 in Merja Zerga (MZ) in the 
Gharb region north of the Sebou River; five in Guercif 
(GCF) on the west side of the Moulouya River and 19 in 
Fritissa (FRT) on the east side of the river in the Oriental 
region (Fig. 1, see Appendix I). We captured the animals 
alive using Sherman traps then we euthanized them by 
cervical dislocation. This protocol was approved by the 
Cuvier (Museum National d’Histoire Naturelle, Paris) 
ethics committee. We extracted a piece of the liver for 
the genetic study and preserved the carcasses in formal- 
dehyde. Later, we extracted and prepared the skulls for 
the morphometric study. Then, we deposited voucher 
specimens in the collections of the Laboratory “Biodi- 
versity, Ecology and Genome’ of the Faculty of Sciences 
of Rabat, except those from the Chrouda locality, that 
are preserved in the collections of the Museum National 
d'Histoire Naturelle of Paris (see Appendix I). 


Morphometric study. We used the standard body mea- 
surements, head-body length (HB), tail length (T), hind 
feet length (HF), ear length (E), weight (WT) and the 
ratio of tail length to head-body length (%T). Using a 
Mitutoyo caliper accurate to 0.01 mm, we took 10 skull 
measurements (mm) on dorsal, ventral and lateral view 
of the skull. We also took three mandibular measure- 
ments (mm). Abbreviations of these craniomandibular 
values are as follows: greatest length of skull (GLS), 
width of the zygomatic arch (WZYG), breadth of brain- 
case (BB), length of nasals (LN), least interorbital con- 


Table 1. Geographic data for the six localities studied in Morocco. 


striction (IO), occipital height from the tympanic bulla 
to the parietal-interparietal suture (HOCC), length of 
anterior palatine foramina (LAF), length of upper molar 
series (M1M3), width of the palate between the first two 
molars M1 at the anterocone (WP), diagonal length of 
tympanic bulla (LTB), length of the mandible from the 
tip of the lower incisor to the angular process (LMDB), 
height of the mandible taken above the mandibular con- 
dyle (HMDB), length of lower molar series (LM1M3). 

The body and skull variability in each population was 
first calculated using descriptive statistics (minimum, 
maximum, mean, standard deviation). Multivariate 
analyses were conducted on log transformed body and 
craniomandibular measurements to make them more ho- 
mogenous. To represent the body and skull variability we 
performed a principal component analysis (PCA). The 
results of the PCA were analysed to identify the factors 
influencing the structure of the data point clouds as visu- 
alized in the bi-plot graphs of the analyses. 

A multivariate analysis of variance (MANOVA) was 
made to test the sexual dimorphism for the body and 
skull variables. We made another MANOVA to test the 
effect of the “locality” variable on the measurements. 
The standard significance threshold for the MANOVA 
was set at 5%. To visualize the morphological variability 
without the influence of sex, we performed a PCA sepa- 
rately within each sex. We also conducted a discriminant 
function analysis (DFA) within each sex to quantify the 
rate of discrimination based on the “locality” variable. In 
this DFA, we did not include the two individuals of CHR 
locality. An additional DFA was conducted on the data of 
the combined sexes to quantify the rate of discrimination 
based on the “locality” variable, and to identify the mea- 
surements that best categorise morphological variabili- 
ty. Finally, we tested the influence of altitude, latitude, 
longitude and annual precipitations using MANOVA 


Locality Altitude (m) Latitude Longitude Precipitations (mm) 
Chrouda 

(CHR) 209 35.360 -5.170 795 
Aounate 

(ANT) 185 32.746 -8.244 317 
Kenitra 

(KTR) 43 34.445 -6.520 579 

Merja Zerga 

(MZ) 25 34.817 -6.310 622 
Guercif 

(GCF) 410 34.206 -3.424 222 
pu 510 34.161 324 218 

(FRT) ; R. 
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(Tab. 1). All statistical analyses were performed using the 
softvvare XLSTAT version 2016 (Addinsoft 2012). 

To minimize age-related bias, two individuals were 
excluded from the statistical analyses because they were 
identified as juveniles. This age estimation was done us- 
ing the dental wear scheme of the upper molar series es- 
tablished by Zyadi (1989) and the age classes based on 


weight made by Zyadi & Benazzou (1992). Also, in all 
the multivariate analyses, we did not include 11 individ- 
uals due to missing data in the measurements. 


Molecular study. We extracted and purified the DNA 
of the 96 individuals using the QIAGEN Kit (DNeasy 
Blood & Tissue Kit) following the manufacturer recom- 


Table 2. Body measurements of Gerbillus campestris in six populations in Morocco. Abbreviations: head-body length (HB); tail 
length (T); hind feet length (HF); ear length (E); weight (WT); ratio of tail length to head-body length (70T), sample size (N); min- 


imum (Min); maximum (Max); standard deviation (S.d.). 


Chrouda Aounate Kenitra Merja Zerga Guercif  Fritissa All popula- 

(CHR) (ANT) (KTR) (MZ) (GCF) (FRT) zu 

N 2 12 24 24 5 19 86 

HB Min 95 83 83 78 81 89 78 
Max 105 103 103 98 116 108 116 

(mm) Mean 100.00 92.33 94.42 90.87 102.40 98.37 94.60 
S.d 7.07 7.02 5.72 4.75 14.54 6.49 727 

N 2 12 22 24 5 18 83 

v Min 109 99 94 93 106 111 93 
Max 119 117 112 117 149 134 149 
(mm) Mean 11400 110.64 106.27 107.54 126.80 121.3 112.29 
S.d 7.07 5.16 5.07 5.59 18.40 7.07 10.09 

N 2 12 24 24 5 19 86 

HF Min 26 23 23 23 24 25 23 
Max 28 26 25 26 29 28 29 

(mm) Mean 2700 24.67 24.52 24.44 25.60 26.55 25.09 
S.d. 1.41 0.78 0.60 0.71 2.07 0.88 1.23 

N 2 12 24 24 5 19 86 

i Min 12 ie 14 13 11 15 11 
Max 15 16 16 16 16 18 18 

(mm) Mean 13.50 14.92 15.10 14.56 14.00 16.18 15.06 
S.d 222 1.16 0.77 0.68 1.87 0.80 1.14 

N 2 12 24 24 5 19 86 

wr Min 34 16.8 17.8 16 18 22 5 16 
Max 40 32 35 25 37.5 36.5 40 

(2) Mean 3700 2397 24.85 20.57 29.00 28.13 24.75 
S.d 4.24 5.37 4.52 2.55 8.31 4.66 5.60 

N 2 12 22 25 5 18 83 

Min 104 111 97 99 H7 104 97 

%T Max 125 131.5 131 144 131 138 144 
Mean 114.54 10165 11:02 118.72 123.96 123.72 119.08 

S.d. 15.17 6.98 9.17 8.77 5.96 10.03 9.97 
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mendations. Then we amplified the cytochrome b gene 
(1040 bp) via polymerase chain reaction (PCR) using 
the primers L7 (ACC AAT GAC ATG AAA AAT CAT 
CGT T) and H15915 (TCT CCA TTT CTG GTT TAC 
AAG AC) (Ducroz et al. 2001). The PCR included an 
initial denaturation step of 3 min at 94°C, followed by 
38 cycles of 30 sec at 94°C, 40 sec at 52°C, and 90 sec at 
72°C, with a final extension step of 5 min at 72°C. Dou- 
ble-stranded PCR products were purified and sequenced 
by Eurofins. We checked the chromatograms then we 
corrected and aligned the sequences both manually and 
using ClustalW in the software BioEdit (Hall 1999). To 
attain the genetic identification, we entered the sequenc- 
es into the Basic Local Alignment Search Tool (BLAST, 
https://blast.ncbi.nlm.nih.gov/Blast.cgi). All the se- 
quences were submitted to GenBank (see Appendix I). 

To choose the model of nucleotide substitution for our 
phylogenetic analysis, we used the software jModeltest 
2.1.10 (Darriba et al. 2012). The General time reversible 
(GTR)+I+G model (Gu et al. 1995) was determined as 
the best-fit model of nucleotide substitution, according 
to the Akaike information criterion (Akaike 1973). We 
constructed a phylogenetic tree using the Maximum 
Likelihood method (ML) in the software MEGA 7.0.26 
(Tamura et al. 2013). The robustness of the obtained to- 
pologies was estimated in all the treatments using a boot- 
strap analysis (1000 replicates). We included in the mo- 
lecular analysis all the cytb sequences of G. campestris 
present in the GenBank database (108 sequences). We 
removed the beginning and the end of all the sequences 
used in the analysis because these parts of the sequences 
were missing in many specimens from GenBank, and we 
maintained a fragment of 1033 pb. We rooted the phy- 
logenetic tree with one representative from each of the 
three subgenera of Gerbillus: G. simoni for the subge- 
nus Dipodillus, G. hesperinus for the subgenus Gerbil- 
lus and G. henleyi for the subgenus Hendecapleura (see 
Appendix I). We also used the Median Joining algorithm 
through the software NETWORK version 4.500 (Bandelt 
et al. 1999) to estimate evolutionary relationships among 
the haplotypes. Finally, the Pairwise Kimura two-param- 
eter (K2P) genetic distances (Kimura 1980) between the 
lineages were computed using MEGA 7.0.26 (Tamura 
et al. 2013). 


RESULTS 


Morphometric study. In the descriptive statistics of the 
body and craniomandibular measurements, we observe 
the existence of a variability between the six populations 
(Tabs 2-3). In general terms concerning body and cranio- 
mandibular size, the highest mean values belong to the 
populations of CHR, GCF and FRT; while the lowest 
mean values belong to the populations of ANT, KTR and 
MZ. 
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Fig. 2. Correlation circle obtained by the principal component 
analysis (PCA) of the body and craniomandibular variables in 
the factor space of (F1 and F2). 


Observations (axes F1 and F2: 74,59 %) 
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Fig. 3. Projection of individuals on the first two axes (F1 and 
F2) of the principal component analysis (PCA). Different co- 
lours represent different localities (Chrouda: black; Aounate: 
green; Kenitra: blue; Merja Zerga: red; Guercif: orange; Fri- 
tissa: purple). 


For the PCA, we visualize in the correlation circle, 
the first two axes F1 and F2 that represent respectively 
64.52% and 10.08% of the total variance (Fig. 2). The 
projection of the 18 body and craniomandibular variables 
on the F1xF2 plane shows that the F1 axis is positively 
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Table 3. Craniomandibular measurements of Gerbillus campestris in six populations in Morocco. Abbreviations: greatest length of 
skull (GLS), width of the zygomatic arch (WZYG), breadth of braincase (BB), length of nasals (LN), least interorbital constriction 
(IO), occipital height from the tympanic bulla to the parietal-interparietal suture (HOCC), length of anterior palatine foramina 
(LAF), length of upper molar series (M1M3), width of the palate between the first two molars M1 at the anterocone (WP), diagonal 
length of tympanic bulla (LTB), length of the mandible from the tip of the lower incisor to the angular process (LMDB), height of 
the mandible taken above the mandibular condyle (HMDB), length of lower molar series (LM1M3), sample size (N); minimum 
(Min); maximum (Max). standard deviation (S.d.). 


Chrouda Aounate Kenitra Merja Zerga Guercif Fritissa All popula- 


(CHR) (ANT) (KTR) (MZ) (GCF) (FRT) tions 
N 2 12 24 24 5 19 86 
Min 28.77 25.42 25.66 25.36 26.98 28.09 25136 
GLS Max 32.66 29 77 28.63 28.98 3107 31.74 32.66 
Mean 30.71 2700 27.43 2721 29.68 30. 12 28.23 
S.d. 705 1.36 0.83 0.80 17: 107 162 
N 2 12 24 24 5 19 86 
Min 12.5 10.77 10.81 10.01 10.97 12 10.01 
WZYG Max 13.8 12.64 12.4 1225 13.26 15:77 13.8 
Mean 13215 11.81 11.62 129 12:32 12.83 11.90 
S.d. 0.92 0.61 0.47 0.49 0.97 0.53 0.81 
N 2 12 24 24 5 19 86 
Min 13.82 12.83 1257 12:59 13.34 13.69 12:57 
BB Max 14.53 14.18 15 15.65 15.59 15.27 15.59 
Mean 14.17 1392 BB 1517 14.57 14.24 1 
S.d. 0.50 0.42 0.26 0.33 0.85 0.41 0.64 
N 2 12 24 24 5 19 86 
Min 11.26 10.09 10.09 9.84 10.43 11.08 9.84 
LN Max 13.51 12.12 1? 17 IL”. 12.87 13.08 1351 
Mean 12.38 1127 11415 11.00 11.89 12.32 11.45 
S.d. 159 0.67 0.49 0.48 1.10 0.49 0:79 
N 2 12 24 24 5 19 86 
Min 5.24 4.52 4.62 4.55 4.88 5.01 4.52 
IO Max 5.8 5.43 52 5.46 5:63 5.94 5.94 
Mean 552 5.00 4.91 4.99 529 5.43 5.10 
S.d. 0.40 0.27 0.14 0.20 0.27 025 0.30 
N 2 12 24 24 5 19 86 
Min 11.2 10.42 9.94 10.05 10771 1111 9.94 
HOCC Max 11.44 11.06 10.87 10.89 1191 11.94 11.94 
Mean 11.32 10.82 10.44 10.44 11.24 11.43 10.78 
S.d. 0.17 0.21 0.26 0.26 0.47 0.23 0.49 
N 2 12 24 24 5 19 86 
Min 52 444 4.2 4.42 4.38 4.9 4.2 
LAF Max 5.64 5.29 5.69 5.65 20 6.08 6.08 
Mean 5.42 4.97 5.23 315 5.16 5.46 522 
S.d. 0.31 029 0.35 0.31 0.48 0:32 0.36 
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Table 3. (continued) 


Chrouda Aounate Kenitra 
(CHR) (ANT) (KTR) 
N 2 12 24 
Min 4.21 3.71 3.38 
M1M3 Max 4.22 4.2 3.89 
Mean 4.21 3.89 3.64 
S.d. 0.01 0.14 0.12 
N 2 12 24 
Min 3.02 2.66 2.88 
VVP Max 3.53 57. 3.54 
Mean 327 2.99 3:19 
S.d. 0.36 0.17 0.16 
N 2 12 24 
Min 8.11 8.14 8.14 
LTB Max JI 8.95 8.91 
Mean 8.60 8.71 8.44 
S.d. 0.70 0.25 0.18 
N > 12 24 
Min 16.37 14.42 14.77 
LMDB Mar 18.56 17.18 16.51 
Mean 17.46 16.05 155: 
S.d. 1.55 0.89 0.50 
N 2 12 24 
Min 6.67 6.31 623 
HMDB Max 7.81 7.49 7.44 
Mean 7.24 6.92 6.85 
S.d. 0.81 0.37 0.28 
N 2 12 24 
Min 3.75 359 327 
LM1M3 Max 4.05 4.03 3.66 
Mean 3.90 3.81 3.53 
S.d. 0:21 0.13 0.10 


correlated with all the variables, which indicates a size 
axis. On the other hand, several variables are positive- 
İy correlated with the F2 axis, while several others are 
negatively correlated with this axis. The most strongly 
positively correlated variable with the F2 axis is the vari- 
able WP, with a correlation value of 0.69. The scatter plot 
of the PCA (Fig. 3) does not allow to visualize different 
distinct eroups, but all the specimens from FRT and CHR 
and virtually all those from GCF are located on the pos- 
itive side of the F1 axis. Contrarily, the specimens from 
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Merja Zerga Guercif Fritissa All popula- 
(MZ) (GCF) (FRT) tions 
24 5 19 86 
3.41 3:97 3.76 3.38 
3.86 4.27 4.36 4.36 
3.65 4.13 4.07 3.81 
0.11 0.11 0.15 0.24 
24 5 19 86 
275 2.85 27 2.66 
3.54 3.41 5 3.58 
3:19 3.18 3.23 3 17 
hir 0.22 0.18 0.19 
24 5 19 86 
7.89 8.55 8.85 7.89 
8.87 9.87 9 75 9.87 
8.44 925 9.36 8.73 
0.24 0.49 0:27 0.46 
24 5 19 86 
14.71 15 16.26 14.42 
16.8 18.35 18.35 18.56 
15:71 1717 1713 16.22 
0.52 1.38 0.63 0.92 
24 5 19 86 
6.42 6.58 7.04 6.23 
7.61 7.74 8.26 8.26 
6.94 7.28 751 7.07 
0.30 0.50 0.36 0.42 
24 5 19 86 
33 3.8 3.74 327 
3.64 4.15 4.13 4.15 
3.50 577 392 3.68 
0.11 0.12 0.10 022 


MZ and KTR are located mainly on the negative side of 
F1. Specimens from ANT show high size variation. 
Among the 83 individuals (55% males, 45% females), 
the average size of females seems somewhat larger than 
the one of males, and the Wilks’ test in the MANOVA 
showed the existence of sexual dimorphism (Wilks’ À = 
0.598, F statistics = 2.386, DF, = 18, DF, = 64, p< 0.05), 
which means that sex has a slightly significant influence 
on the dispersal of the observations on the PCA. For the 
PCA of the male individuals, we visualize the scatter plot 
of the two first axes F1 and F2 that represent respective- 
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Fig. 4. Correlation circle obtained by the discriminant function 
analysis (DFA) of the body and craniomandibular variables in 
the factor space of (F1 and F2). 


ly 66.45% and 10.17% of the total variance (see Appen- 
dix IT). For the PCA of the female individuals, we visu- 
alize the scatter plot of the two first axes F1 and F2 that 
represent respectively 62.55% and 13.46% of the total 
variance (see Appendix III). In both the PCA of males 
and females, the F1 axis is positively correlated with 
all the variables indicating a size axis; with the cranio- 
mandibular measurements HOCC, LTB, WZYG, BB, 
LN, HMDB, LMDB and GLS being the most strongly 
correlated variables with F1 (values ranging from 0.83 
to 0.98). In both scatter plots (males and females), indi- 
viduals from MZ and KTR are the smallest, while those 
from FRT are the largest. Male individuals from GCF are 
almost as large as the individuals from FRT, while the 
female individuals from GCF have a high size variation. 
Male individuals from ANT show a medium size and the 
females show high size variation. The male individual 
from CHR 1s one of the largest individuals and the female 
from CHR has a medium size. In the DFA, we were able 
to quantify the discrimination between five populations 
within each sex: in the males, 100% of the individuals of 
ANT were correctly classified, 86.67% for KTR, 81.25% 
for MZ, 100% for GCF and 100% for FRT; while in the 
females, 100% of the individuals were correctly classi- 
fied in each one of the five populations. 

The Wilks’ test in the MANOVA revealed the presence 
of a highly significant morphological difference between 
the various localities (Wilks’ À = 0.01, F statistics = 
5.146, DF, = 90, DF, = 296, p < 0.0001). According to 
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this result, we performed a discriminant function analysis 
(DFA) to quantify the morphological variability based on 
a locality effect. In the correlation circle (Fig. 4), the first 
two axes Fl and F2 of the DFA represent respectively 
70.68% and 10.68% of the total variance. The projec- 
tion of the 18 body and craniomandibular variables on 
the F1xF2 plane shows that the F1 axis is positively cor- 
related with all the variables except WP (-0.02). In the 
F1 axis, the most discriminating variables are LM1M3 
(0.91), HOCC (0.88), M1M3 (0.87), LTB (0.82), BB 
(0.80) and GLS (0.75). 


Observations (axes F1 and F2: 81,36 76) 


F2 (10,68 96) 


_F1 (70,68 %) 


Fig. 5. Projection of individuals on the first two axes (F1 and 
F2) of the discriminant function analysis (DFA). Different sym- 
bols represent different localities. Also shown are confidence 
circles for each locality using covariance hypothesis. 


On the scatter plot of the two first axes of the DFA 
(Fig. 5), although some overlap can be observed, we can 
distinguish four major groups: 


— Group I composed of the individuals of Merja Zerga 
(MZ) and Kenitra (KTR); 

— Group II composed of the individuals of Aounate 
(ANT); 

— Group III composed of the individuals of Guercif 
(GCF) and Fritissa (FRT); 

— Group IV composed of the individuals of Chrouda 
(CHR). 


This structuration is integrated by the F1 axis that iso- 
lates Group I on the negative side of the axis and Groups 
III and IV on the positive side, while Group II is located 
in an intermediate position along this axis. The F2 axis 
allows a complete separation between Group IV and 
Group III. The DFA has also allowed to quantify the dis- 
crimination between the six populations: 100% of the 
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Fig. 6. Phylogenetic tree of Gerbillus campestris for the cytochrome b gene resulting from the Maximum-Likelihood analysis 
(GTR + I + G substitution model). Numbers at nodes represent ML bootstrap support. To improve clarity, values of the most apical 
nodes are not shown. The scale bar represents the branch length measured in the number of substitutions per site. The ten mitochon- 
drial lineages are presented on the right border along with their geographic origin. 
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individuals of ANT were correctly classified, 86.4% for 
KTR, 87.5% for MZ, 80% for GCF, 100% for FRT and 
100% for CHR. 

Based on the results of the MANOVA, all the geo- 
graphic parameters (altitude, latitude, longitude, annual 
precipitations; Tab. 1) have a highly significant effect: 
Wilks’ A = 0.01, F statistics = 5.146, DF, = 90, DF, = 
296, p < 0.0001). A supplementary observation concern- 
ing the body variation between our suggested morpho- 
metric groups is the difference of dorsal coat colour, 
being dark brown for Group 1, II and TV, and yellowish 
brown for Group III. 


Molecular study. Based on the results of the phyloge- 
netic tree (Fig. 6) and the haplotypes network (Fig. 7), 
the individuals from the six localities belong to four of 
the nine phylogenetic lineages previously identified by 
Nicolas et al. (2014): 


— All the ten individuals from Chrouda (CHR) are found 
in lineage 7, which was previously composed of four 
specimens from this same locality. 

— All the 13 individuals from Aounate (ANT) are found 
in lineage 1, which is widely distributed in Morocco 
and contains specimens from many localities south of 
the Bou Regreg River. 

— All the 24 individuals from Kenitra (KTR) and the 25 
individuals from Merja Zerga (MZ) are found in lin- 
eage 6, that was previously composed of only nine 
specimens from Merja Zerga. 

— All the five individuals from Guercif (GCF), west of 
Moulouya River, and the 19 individuals from Fritissa 
(FRT), east of Moulouya River, are found in lineage 4 
that previously contained only specimens from the east 
of this river (Morocco, Algeria and Tunisia). 


Moreover, we have added three individuals from Egypt 
(Khalifa et al. 2018) and one from Libya (Chevret et al. 
2014) to lineage 2, that was previously composed only 
of specimens from Niger and Mali (Tessalit and Adrar 
des Iforas localities). We also added one individual from 
Tunisia (Alhajeri et al. 2015) to lineage 4. Furthermore, 
we have attributed one individual from Mali (Sinkerma 
locality) (Schwan et al. 2012) to a new suggested lineage 
10. The K2P genetic distance between lineage 10 and the 
other lineages goes from 1.4% with lineage 8, to 2.6% 
with lineage 2. 


DISCUSSION 


Morphological adaptation. The results of the morpho- 
metric study revealed a significant morphological vari- 
ability among Gerbillus campestris in Morocco. The 
significant difference in the MANOVA testing the effect 
of “locality” variable; and the relatively high rate of dis- 
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crimination of the DFA indicates a spatial structuration 
between populations. This structuration was represented 
by Group I (Merja Zerga and Kenitra populations), Group 
II (Aounate population), Group HI (Guercif and Fritissa 
populations) and Group IV (Chrouda population). The 
parameters that contribute the most to this clustering are 
the lower and upper molar series length, the skull height, 
the tympanic bulla length and the braincase breadth. The 
variation of these parameters shows a gradual increase in 
size, from Group 1 being the smallest, to Group TV and 
Group III being the largest; while Group II has a medium 
size. In Morocco, a few previous studies were made on 
the geographic variation in biometrics for G. campestris 
(Benazzou & Zyadi 1990; Baala 1995). Apart from the 
Kenitra locality, they showed for different localities than 
ours an interpopulational variation in size, with Kenitra 
being the smallest. Another type of body variation re- 
ported in Morocco was related to the dorsal coat colour, 
being fawn grey in the northern regions and gradually 
brightening towards the south, until becoming light fawn 
in desertic habitats (Petter 1961). 

In our study, the morphological variability was slightly 
influenced by sex, as we showed in the MANOVA the 
existence of a weak sexual dimorphism when males and 
females were analysed together. The PCA of the variabil- 
ity within each sex shows almost similar results between 
males and females; and roughly a similar distribution of 
individuals in the F1xF2 plane. Also, there is a strong rate 
of discrimination between the populations within each 
sex. All these reasons support a weak influence of sex 
on the clustering of our populations in distinct groups. In 
the precedent study of G. campestris in Morocco, Baala 
(1995) found no significant sexual dimorphism in five 
populations. In Tunisia, Jordan et al. (1974) found sex- 
ual difference in only one character. Additionally, sexual 
dimorphism has seldom been confirmed in Muridae ro- 
dents using classical morphometric analysis (Csanady & 
Mo$ansky 2018). Thus, larger samples from each pop- 
ulation are needed to verify if the minor sexual dimor- 
phism that we uncovered is not due to sampling bias. 

The morphometric structuration that we showcased in 
G. campestris is explained by the wide distribution of 
this species in Morocco and the vast environmental di- 
versity of this country. The significant differences found 
in the MANOVA for the geographical factors (altitude, 
latitude, longitude, precipitations) give more explanation 
to support this spatial and environmental differentiation. 
A similar result was found for Rattus rattus in Tunisia, 
where latitude, longitude, altitude and precipitation vari- 
ables exhibited a significant control on size parameters 
(Ben Faleh et al. 201 2a). 

Our different populations are located in diverse bio- 
climates: subhumid for Group I (Kenitra and Merja 
Zerga) and Group IV (Chrouda), semi-arid for Group II 
(Aounate) and arid for Group III (Guercif and Fritissa) 
(Mokhtari et al. 2013). The influence of varying biocli- 
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Fig. 7. Minimum spanning network of Gerbillus campestris for the cytochrome b gene haplotypes. Circle sizes are proportional 
to the number of similar haplotypes observed in the data set. Branch lengths are proportional to the number of mutations between 
haplotypes. The dashed ellipses (added for clarity) represent the ten mitochondrial lineages. 
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mates and rainfall patterns on the populations acts indi- 
rectly through the availability and type of food resources 
(Yom-Tov öz Geffen 2011, Breno et al. 2011). Previous 
studies shovved an influence of climate on the skull of 
rodents, through variations in the type and availability of 
food (Renaud et al. 2005; Samuels 2009). The quality and 
quantity of food were shown to influence growth rates 
and skull morphology in the multimammate rat, Mas- 
tomys natalensis (Lalis et al. 2009; Breno et al. 2011). In 
the muskrat, Ondatra zibethicus, the high availability of 
food to each individual accelerated its growth, but if food 
was less abundant, there was an extension of the growth 
phase and an increase in growth variability (Pankakos- 
ki et al. 1987). Pocket gophers inhabiting alfalfa crops 
had larger body and skull size compared to those living 
in natural dry environments (Patton & Brylski 1987). In 
another example, a high difference in protein content in 
the diet of captive rats, Rattus norvegicus, caused a clear 
modification in adult skull shape but much less modifi- 
cation in the adult body and skull size (Miller & German 
1999; Reichling & German 2000). However, squirrel 
monkeys, Saimiri sciureus that were fed a diet varying in 
protein content had different skull sizes but similar shape 
(Ramirez Rozzi et al. 2005). Surprisingly in our study, 
the individuals of Kenitra and Merja Zerga are the small- 
est even if they occupy the fertile agricultural fields of 
the Gharb region, characterized by high food availability 
and the abundance of peanut crops known to be rich in 
protein; as opposed to the individuals of Guercif and Fri- 
tissa, that are the largest even if they occupy habitats with 
low food availability in the Oriental region. 

Food hardness, too, has an influence on skull devel- 
opment. This is an important factor of differentiation 
between our suggested morphometric groups, especially 
that the most discriminating variables in the DFA are the 
lower and upper molar series length. It was established in 
various species of rodents that skull morphology can be 
associated with the structure of food. In the prairie deer 
mouse, Peromyscus maniculatus, a softer diet caused a 
reduction in size of the masseteric tubercles and a nar- 
rowing of the zygomatic plate (Myers 1996). In laborato- 
ry rats, individuals feeding on a soft diet suffered a 12% 
loss of the mass of the skull and mandible, 4% decrease in 
the length of the angular process and 1% to 2% decrease 
in the dimensions of the cranial and facial skeleton, com- 
pared to the individuals feeding on a normal diet (Moore 
1965). Other studies performed on rats demonstrated that 
maxillary breadth was markedly increased in the hard-di- 
et reared individuals (Beecher & Corruccini 1981); that a 
modification in food consistency altered the masticatory 
muscle function and thus had an impact on craniofacial 
morphology (Kiliaridis 2006); and that a harder diet in- 
creased the growth of craniomandibular structures, espe- 
cially those connected to the masticatory function (Abed 
et al. 2007). Feeding on a harder diet also relatively aug- 
mented the skull width of other mammals (Corruccini & 
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Beecher 1984; Ciochon et al. 1997). In our studied pop- 
ulations, the smallest individuals, representing Group I 
populations (Kenitra and Merja Zerga) feed on various 
crops but primarily on peanut crops, in the form of seeds 
at the seedling stage, or the plant at the vegetative stage, 
or the pods at maturity (Zyadi & Benazzou 1992). On the 
opposite, the biggest individuals, representing Group III 
(Guercif and Fritissa), feed mainly on wheat and other 
cereals, in the form of plants and their hard grains. Our 
explanations connecting the morphological differences 
to the quantity and quality of food resources are hypo- 
thetical but reasonable, even if we do not possess quanti- 
fiable data on the food resources available to the studied 
populations or the stomach content of the individuals. 

Differentiation between our morphometric groups was 
also partly due to the increase in size of the occipital part 
of the skull, represented by an increase in the size of the 
tympanic bulla, and therefore an increase in skull height 
and braincase breadth. This hypertrophy is observed in 
our populations with the increase of aridity, from Group I 
with the smallest bullae, found in a subhumid bioclimate, 
to Group III with the largest, found in an arid bioclimate. 
The size increase of the bullae also seems to go along 
with the decline of the population’s density, from Group 
I of the Gharb region where G. campestris is prolific 
and can pullulate (Zyadi & Benazzou 1992), to Group 
III of the Oriental region where this species is much less 
abundant. Bioclimate variability, through the change of 
aridity levels, causes variations in vegetation cover and 
the abundance of food resources. At first, fewer food re- 
sources are available in the more arid habitats, and con- 
sequently there is a lowering of the rodent population 
density. Secondly, more arid habitats are characterized 
by a lower vegetation cover and thus more open habitats. 
Hypertrophy of the tympanic bullae, and consequently 
enhanced auditory sensitivity, Is considered as an adap- 
tation to these conditions of arid habitats. Gerbillinae ro- 
dents were distinguished to be one of the rodent families 
most predisposed to the hypertrophy of their bullae (Pet- 
ter 1961). A better hearing acuity facilitates the encounter 
of the sexes and reproduction in low population densities 
(Petter 1961) and can help in evading predators in these 
more open habitats (Lay 1972; Webster & Webster 1975; 
Alhajeri et al. 2015). Similar to our results, Colangelo 
et al. (2010) found in Gerbilliscus that modifications in 
the molar row and the tympanic bulla were related to cli- 
mate. 

Other possible factors influencing morphological vari- 
ability are predation pressure and intra- and interspecific 
competition (Yom-Tov & Geffen 2011). The different 
degrees of predation pressure by various predators (birds 
of prey, reptiles, small carnivores, rats, dogs) as well as 
competition over resources with the other rodents that we 
have captured in the studied localities (Mus musculus, 
Mus spretus, Rattus rattus, Rattus norvegicus, Meriones 
shawii), need to be tested by joint ecological studies. 
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Phylogeographic distribution. Previously, Nicolas et al. 
(2014) revealed a complex phylogeographic structuration 
of Gerbillus campestris in the form of nine mitochondrial 
lineages in Morocco, North Africa (Algeria and Tunisia) 
and the western Sahel region (Mali and Niger). Accord- 
ing to these authors, different factors contributed to this 
structuration: isolation by distance; the climatic fluctua- 
tions of the Pleistocene marked by warming and cooling 
conditions, and the period of aridification at the end of 
the Pleistocene; and the role of the Moulouya River and 
Bou Regreg River as biogeographical barriers to genetic 
flow. The new sampled specimens and populations in our 
study give additional details concerning the distribution 
of the various lineages (Fig. 1). 

The supplementary individuals from Chrouda locali- 
ty reinforce the presence of lineage 7 in this locality as 
an allopatric lineage, in spite of the small geographic 
distance between this locality and the Esperada locality 
where exists the allopatric lineage 8, that forms a recent 
monophyletic group with lineage 7 beyond the Rif moun- 
tains in the extreme north of Morocco. 

The added individuals from Aounate locality show that 
lineage 1 is still the only lineage present in this locality, 
even though this lineage is sympatric with lineage 3 in 
almost all the localities surrounding Aounate locality, in 
the north, east and south; between the Bou Regreg River 
and Souss River. 

In the allopatric lineage 6, the Merja Zerga locality and 
the new locality of Kenitra are both located between the 
northern side of the Sebou River, the largest river in Mo- 
rocco, and the southern side of the Loukkos River. How- 
ever, on the southern side of the Sebou River, in the Sidi 
Boughaba locality just nearly 30 km from the Kenitra lo- 
cality, we find the allopatric lineage 5. This suggests that 
the Sebou River is a significant barrier to gene flow for 
this species. Studies on different plant species have found 
that the Sebou River was a major geographical limit to 
gene flow (Hypochaeris salzmanniana DC. by Ortiz 
et al. 2007; Chamaerops humilis L. by Garcia-Castafio 
et al. 2014; Onopordum dissectum Murb. by Balao et al. 
2017). The north-west of Morocco is surrounded by the 
Rif and Middle Atlas Mountains and comprises several 
long rivers (Bou Regreg, Sebou and Loukkos Rivers) 
originating in these mountains and extending to the At- 
lantic coast at different latitudes, creating large estuaries 
and floodplains (Chichagov 2008), that probably acted as 
barriers separating the populations of G. campestris and 
promoting genetic isolation. 

Lineage 4 was thought to be restricted to the eastern 
bank of the Moulouya River in Morocco and extend- 
ing to the north of Algeria and to the north of Tunisia 
(Nicolas et al. 2014). However, thanks to our sampling 
in the new localities of Guercif and Fritissa, we showed 
that lineage 4 originating on the eastern bank of the 
Moulouya has later spread towards the western bank of 
the river. For another North African Gerbillinae rodent, 
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Meriones shawii, Lalis et al. (2016), using both mtDNA 
and nDNA, suggested to split it into two species with 
an east-west structuration roughly around the Moulouya, 
though samples were missing from the western side of 
the Moulouya and from Northern Morocco. On the con- 
trary, in the Greater Egyptian Jerboa, Jaculus orientalis, 
one of the three identified cytb lineages is found both in 
northern Morocco west of the Moulouya and in north- 
ern Algeria (Ben Faleh et al. 2012b). Based on Beddek 
et al. (2018), the structuration of mitochondrial diversity 
of reptiles and amphibians species around the Moulouya, 
between northern Morocco and north-western Algeria, 
can be assembled into three groups: 1) the same lineage 
exists in north-western Algeria and northern Morocco, 
on both sides of the Moulouya (Timon tangitanus, Tro- 
gonophis wiegmanni, Podarcis vaucheri, Acanthodacty- 
lus erythrurus); ii) north-western Algeria is inhabited by 
the same lineage that is found in northern Morocco but 
only on the eastern side the Moulouya (Natrix maura), 
111) a continuity between Moroccan populations from the 
two sides of the Moulouya, but with a separation from 
the populations of north-western Algerian (Pelophylax 
saharicus, Natrix maura). For these reptiles and amphib- 
ians, the origin of the east-west phylogeographic break 
around the Moulouya was attributed to the former marine 
Rifian corridor that joined the Atlantic Ocean with the 
Mediterranean Sea, east of the Rif mountains, between 
8.0 Mya and 5.6 Mya (Barhoun & Bachiri 2008; Beddek 
et al. 2018). Although part of this marine corridor was 
located on the present Moulouya valley, this paleogeo- 
graphic barrier cannot be applied, at least not directly, to 
the case of G. campestris since the history of this species 
has occurred starting from the middle Pleistocene (0.4- 
0.5 Mya) (Tong 1989). Our results may lessen the role of 
the Moulouya River and the valley encompassing it, in 
preventing the spread of G. campestris towards western 
Morocco at some point in the past, but the real geograph- 
ic differentiation in this region 1s complex and should be 
further explored. An alternative possible explanation for 
the occurrence of lineage 4 on the western bank of the 
Moulouya River is a human-mediated transport of indi- 
viduals, which can represent a counterbalancing factor to 
the genetic differentiation caused by geographic barriers. 

Surprisingly, and although they originate from very 
distant localities, the specimens from Libya and Egypt 
belong to the Saharan lineage 2 along with the individ- 
uals of Niger and Mali (Tessalit and Adrar des Iforas 
localities), which makes lineage 2 the most widespread 
lineage, being present in Mali and Niger in the south and 
west; and Libya and Egypt in the north and east. Ger- 
billus campestris has a large distribution and occupies 
a wide variety of habitats, except in the Sahara Desert 
where it has a sporadic distribution and is restricted to 
rocky massifs. The Sahara Desert acts as a barrier to 
genetic flow between north and south of the desert, for 
species that are not adapted to extremely arid and san- 
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dy habitats, like G. campestris. However, the presence 
of lineage 2 from Mali and Niger in the south, to Egypt 
and Libya in the north of the Desert can represent a relic 
of a larger distribution during more favourable periods 
of the Pleistocene. Paleoenvironmental data and climatic 
reconstructions show that North Africa regularly shifted 
from wetter to drier climatic conditions throughout the 
Quaternary (deMenocal 2004; Schuster et al. 2006; Ta- 
bel et al. 2016), leading to periodic contraction/expan- 
sion of the Sahara Desert. Similarly, Nicolas et al. (2018) 
showed through ecological niche modelling and genetic 
data that presently unsuitable arid areas of the Saharan 
Desert were suitable for the amphibian Bufotes bouleng- 
eri throughout most of the Quaternary. 

Another unexpected result concerns the specimen from 
south central Mali (Sinkerma locality) representing the 
added lineage 10. Even if it is geographically close to the 
specimens from north-east Mali (Tessalit and Adrar des 
Iforas localities), the Sinkerma individual is genetically 
very distant from lineage 2 (2.6%) and seems to be ge- 
netically closer to lineages from the north of Morocco. In 
this case, the high genetic differentiation between the ad- 
jacent populations eliminates the role of isolation-by-dis- 
tance as a cause of genetic variability. The Sinkerma 
locality is found south of the Niger River, the third lon- 
gest river in Africa, compared to the two other localities 
that are situated north of the river. Therefore, the Niger 
River may have played a role as a barrier to gene flow 
between lineages, similarly to what was suggested for 
other rodents (Nicolas et al. 2008, 2009; Brouat et al. 
2009; Dobigny et al. 2010; Hima et al. 2011; Colangelo 
et al. 2013). Also, phylogeographic studies for different 
rodents showed the existence of distinct phylogroups in 
this region of West Sahel (Mastomys erythroleucus by 
Brouat et al. 2009; Acomys chudeaui by Nicolas et al. 
2009; Gerbillus henleyi by Bouarakia et al. 2018). To 
explain the phylogeographic discontinuity for G. camp- 
estris in Mali, more samples from this region are needed 
to evaluate the combined actions of barriers and allopat- 
ric differentiation or the effects of lineage sorting from 
a highly polymorphic ancestral gene pool (Avise 2000). 


Relationship between morphometric and molecular 
variability. A crucial goal of this work is the confronta- 
tion of the results of the morphometric study to those of 
the molecular study. Previously, Baala (1995) showed in 
five populations in Morocco, that the biometric variabil- 
ity was not accompanied by allozymic variability based 
on biochemical differentiation of proteins. However, in 
our study, the morphometric results are congruent with 
the genetic results and the geographic distribution. The 
geographic structuring of the six populations of G. camp- 
estris in Morocco has produced sufficient phenotypic dif- 
ferentiation and high levels of sequence divergence that 
allows their discrimination into four groups: 
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— Group I of the morphometric study, composed of the 
populations of Merja Zerga and Kenitra, corresponds 
to the genetic lineage 6 and is present in north-west- 
ern Morocco. 

— Group II, representing the population of Aounate, 
corresponds to the genetic lineage | and is present in 
central Morocco; 

— Group III, composed of the populations of Guercif 
and Fritissa, corresponds to genetic lineage 4 and is 
present in north-eastern Morocco; 

— Group IV representing the population of Chrouda, 
corresponds to the genetic lineage 7 and is present in 
the extreme north of Morocco. 


A notable point is the clarification of the phylogeograph- 
ic status of the populations of Guercif and Fritissa in 
relation to the populations present either on the eastern 
bank or the western bank of Moulouya River. This river 
was considered to be a major barrier reducing dispersal 
rates and thus hindering gene flow between the eastern 
populations (Tunisia, Algeria and east of Morocco) and 
the western ones (west of Morocco). Yet, both the mor- 
phological and genetic analyses showed that the popula- 
tions of Guercif (west of Moulouya) and Fritissa (east of 
Moulouya) belong to the same group. 

Although we did not determine the specific relationship 
between morphological and genetic variation, we were 
able nonetheless to demonstrate for G. campestris its high 
plasticity, its capacity of adaptation and the micro-evolu- 
tionary processes at work. In summary, G. campestris 1s a 
widely distributed species that can be divided in Moroc- 
co into different assemblages of populations that experi- 
enced variable degrees of geographic isolation and adapt- 
ed to different local environmental conditions, resulting 
in a strong morphological differentiation and genetic 
structuration. But systematically speaking, the question 
to be asked is whether we should see this variability as 
a progressive gradient or must divide the species into 
different diverging subspecies. Several subspecies were 
previously described based on geographic variations of 
pelage coloration, external size and/or skull characters 
(Setzer 1958, Petter & Saint-Girons 1965; Ranck 1968, 
Jordan et al. 1974; Osborn & Helmy 1980). Based only 
on their geographic distribution, we may attribute some 
populations used in the current study to these different 
subspecies. For example, the individuals from Tunisia 
and Algeria belonging to lineage 4 may be attributed 
to the subspecies G. campestris campestris Levaillant, 
1857. However, in other cases, specimens from the same 
lineage fall within different subspecies: the individuals 
from Egypt fit with G. campestris haymani Setzer, 1958; 
while the one from Libya fits with G. campestris dod- 
soni Thomas, 1902. Furthermore, Petter & Saint-Girons 
(1965) suggested that only two subspecies are found in 
Morocco and they were identified based merely on color- 
ation criteria (G. c. campestris and G. c. dodsoni), which 
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does not correspond to the complex phylogeographic 
structuration revealed in this country. Gerbillus hilda 
Thomas, 1918, first described as a new species in Mo- 
rocco on the northern Atlantic coast (122 km south-west 
of Tangiers), approached the dimensions of G. campes- 
tris externally and cranially but was characterized by a 
small band of hairs on the hind feet soles (Thomas 1918). 
Specimens captured 17 km south-west of Rabat and de- 
scribed as G. hilda by Schlitter & Setzer (1972) appear to 
be closer in body and skull dimensions to the G. camp- 
estris individuals of Group I. G. hilda was later regarded 
as a subspecies of G. campestris (Petter 1975 and Corbet 
1978, as cited in Lay 1983), while Lay (1983) viewed 
it as a subspecies of G. nanus. If we want to reconsider 
the taxonomy of Gerbillus campestris on an intraspecific 
level, genetic variability and morphological character- 
istics must be further identified across the range of the 
species. 


CONCLUSIONS 


Supplementary sampling from the whole distribution 
area of Gerbillus campestris should be carried out, and 
different markers must be used in the genetic study to 
better explore the genetic diversity of this species in the 
Sahara Desert (Mali, Niger) and in Algeria. In addition, 
the study of morphological variability should be refined 
in order to include specimens from all the geographic 
range of the species and all the genetic lineages. Clas- 
sical morphometry is not sufficient and craniomandibu- 
lar variation should be further studied using geometric 
morphometric methods. Further research should focus on 
how morphological and genetic variations could be relat- 
ed to the adaptation capacity of this rodent to anthropized 
habitats, increasing its pullulation capacity and the dam- 
age it can cause to agricultural crops. 
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APPENDIX IL 


Projection of male individuals on the first two axes (F1 and F2) of the principal component analysis (PCA). Different colours 
represent different localities (Chrouda: black; Aounate: green; Kenitra: blue; Merja Zerga: red; Guercif: orange; Fritissa: pur- 
ple). The dashed ellipses (added for clarity) represent the six populations. 


Observations (axes F1 and F2: 76,63 %) 


F1 (66,45 %) 
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APPENDIX III. 


Projection of female individuals on the first two axes (F1 and F2) of the principal component analysis (PCA). Different colours 
represent different localities (Chrouda: black; Aounate: green; Kenitra: blue; Merja Zerga: red; Guercif: orange; Fritissa: purple). 
The dashed ellipses (added for clarity) represent the six populations. 
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